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Abstract 

This work presents exchange potentials for specific orbitals calculated by invert¬ 
ing Hartree-Fock wavefunctions. This was achieved by using a Depurated Inversion 
Method. The basic idea of the method relies upon the substitution of Hartree-Fock 
orbitals and eigenvalues into the Kohn-Sham equation. Through inversion, the corre¬ 
sponding effective potentials were obtained. Further treatment of the inverted potential 
should be carried on. The depuration is a careful optimization which eliminates the 
poles and also ensures the fullfilment of the appropriate boundary conditions. The 
procedure developed here is not restricted to the ground state or to a nodeless orbital 
and is applicable to all kinds of atoms. As an example, exchange potentials for noble 
gases and term-dependent orbitals of the lower configuration of Nitrogen are calcu¬ 
lated. The method allows to reproduce the input energies and wavefunctions with a 
remarkable degree of accuracy. 
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1 INTRODUCTION 


The successful idea of replacing a many-body, non-local interaction by an effective one- 
electron equation opened up the possibility of studying extremely complex systems with 
high accuracy. In this context, the success of the Kohn-Sham density-functional theory-^^ 
(DFT) began when crucial developments in its exchange-correlation terms gave the the¬ 
ory predictive power to compete with well-developed wavefunction methods^. The impor¬ 
tance of the exchange-correlation potentials in chemical physics has been emphasized by 
Bartlett^'^. Exchange potentials are in general constructed by local approximations to the 
nonlocal Hartree-Fock exchange operator (i.e. the Slater potential-, the optimized effective 
potential-^-, the Krieger-Li-Iafratel^ and several other&ii^— ). 

The atomic collision community, on the other hand, is also eager to accurately determine 
effective one-electron local potentials which would allow to generate in a simpler way the 
wavefunctions of the particles interacting in a scattering process. In particular, we need 
to represent an orthonormal set of bound and continuum states to calculate the transition 
probabilities. This should include detailed n/-orbital potentials, a feature missing in most 
of the standard density functional methods. Soft pseudopotentials like abinit— or USPP— 
cannot be used because they overlook the information of the internal region of the wavefunc¬ 
tions. The features of this region can play a very important role, such as the cusp conditions 
in the processes of electron capture and ionization. In an attempt to meet the needs of 
both chemist and collisionist communities, we strove to obtain accurate and simple specihc 
nZ-orbital local potentials. 

How to determine central potentials from known electron wavefunctions and densities is 
a well studied subject in the DFT community— The extraction of the true Kohn-Sham 
exchange-correlation potential from near-exact electronic densities has been demonstrated, 
with particular reference to two-electron systems like He^, He-isoelectronic ion&^i, and 
H 2 — >2^ as well as exact soluble models (for example, an external harmonic potential as in 
Filippi et al—). 

Some other works start with a particular Kohn-Sham potential and solve the correspond¬ 
ing equations, obtaining the KS orbitals^i^^. Through inversion, they obtain a reconstructed 
KS potential, which agrees almost everywhere with the original one, except in some regions 
where huge oscillations arise. In some cases, the reconstructed potential may be distorted 
beyond recognition— The same type of procedure was suggested many years ago by 
Hilton et al., in applications circumscribed to the calculation of photoionization processes 
of atoms2^>22^ water— and other molecules^*^^. These papers, in turn, refer to the earlier 
work in atomic polarizability carried out by Sternheimer— and Dalgarno and Parkinson—. 
However, they focused on the hnal photoionization cross section results, and did not provide 
details about the quality of the potentials and the wavefunctions they generated. 

Assuming the validity of the separation between exchange and correlation functionals, we 
will focus here only on the calculation of the exchange contribution to the potential. Since 
Hartree-Fock does not include the correlations, our approach allows to obtain the “exact” 
one-electron local potential representing the exchange interactions. Strictly speaking, the 
method does not rely on the KS inversion formula since the Hartree-Fock solutions were the 
ones used for the inversion. That is, we solved a KS-type equation, but rather than having 
KS-orbitals, we operated directly with the Hartree-Fock wavefunctions. For open-shell 
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atoms, we were able to find orbital spin-polarized exchange potentials, this being crucial, 
for instance, to find the hyperfine coupling constants^>^. 

However, this is not a simple task, and probably that is why the method presented here 
has not been widely applied in the past. If the wavefunction has nodes, it will produce huge 
poles in the potential. Moreover, even for nodeless states, the asymptotic decaying behavior 
of the bound wavefunctions produces severe numerical difficulties, making the inversion op¬ 
eration intractable sometimes. In our method, a depuration procedure follows the inversion. 
This depuration implies, first, the annihilation of the poles. Then, a careful optimization of 
the potential which ensures the fulfillment of the appropriate boundary conditions. 

The work is organized as follows. Section |2] describes the method, which includes the 
inversion procedure fl2.ip . the potential depuration fl2.2p and its further optimization fl2.3p . 
Section |3] presents the resulting effective potentials for the orbitals corresponding to the 
ground states of different noble gases, including a thorough examination of the wavefunctions 
generated by these potentials fl3.ip . The corresponding exchange potentials are discussed in 
fl3.2p . comparing the potentials for specific-n/ orbitals with averaged potentials. Results of 
the same calculations for the Nitrogen atom are provided in fl3.3p . Atomic units are used 
unless otherwise specified. 


2 THEORY 


2.1 The direct inversion method 


The radial part of the Schrodinger equation for an electron in a local and central potential is 
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We assume the following hypothesis: If the wavefunctions Uni are replaced by the solutions 
of an Hartree-Fock calculation u^f, then, the corresponding effective local potentials 
that generate such wavefunctions should exist. Based on this we converted the HF method 
into a set of Kohn-Sham equations, whose solutions are the Hartree-Fock wavefunctions: 
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The effective potentials given by. 


VS^{r) = V^{t) + , 


(3) 


are composed of the external potential V'" (the Coulomb field of the nucleus), the direct 
(or Hartree) potential (the electrostatic electron repulsion), and the orbital exchange 
potentials V^i- We have ignored the correlation term since the HF solutions do not include it. 

Since the solutions are known (calculated numerically with the HF code by C. F. 
Fischer—, and the nrhf code by W. Johnson—) we proceeded to directly invert the Kohn- 
Sham-type equations: 



^ ^ d? HF/ N _ + 1) , HF 
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obtaining the HF inverted potential Assnming a Coulombic-type shape, it is con¬ 

venient to dehne an HF inverted charge 

= (5) 

The direct computation of dl]) is known to pose serious numerical problems^. First, the 
presence of (genuine) nodes in the wave function to be inverted produces poles and unrealistic 
features around them. This has led to the general consensus that the inversion method can 
only be used for nodeless orbitals^. Second, numerical rounding up of the exponential decay 
of the bound states hinders the corresponding inverted potential from having the physically 
desired asymptotic form. Moreover, there is a third problem at the very heart of the Hartree 
Fock method: the exact solutions may have oscillations (and therefore, spurious nodes) in 
the large-r or “tail” region of the functions. The existence of these spurious nodes in Hartree 
Fock was already suggested by Fischer—. This failure is not caused by the numerical scheme 
but it is inherent to the method. Probably, these nodes are surviving long-range exchange 
effects due to the non-local character of the Hartree-Fock wavefunctions: the behavior of a 
particular orbital depends on all others. We have found the same spurious nodes at the same 
places even using different numerical codes. As a general rule, the spurious nodes appear 
at very long distances, in regions where the amplitude of the wavefunction is very small. 
Therefore, their existence has no practical consequences, and they can be ignored in any 
general Hartree-Fock calculation. However, this is not true as far as the inversion procedure 
is concerned, as we will discuss in the next section. Other examples where the presence 
of orbital nodes (both formal and those in the tail region) can be problematic in inversion 
procedures can be found in the literature (see for instance Peach et al.—). 


2.2 The depurated inversion method 

The difficulties mentioned above make it very hard to obtain the correct V^^{r) potentials 
using the simple inversion formula given by Eq. (0]). To overcome these troubles we have 
developed a depurated inversion method (DIM) which optimizes the effective charges rather 
than the effective potentials. We managed to constrain any potential to have the right 
boundary conditions by enforcing the effective depurated inverted charge to behave as follows: 


ryDlM 

^nl 
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Zjv as r —)■ 0 

1 as r —)■ cxD 


( 6 ) 


where Zn is the nuclear charge. Once the charge is determined at the boundaries, we can 
obtain a smooth analytic expression for htting the Z^^ (r) for the largest possible 

range, except in the neighborhood of the nodes. All this can be accomplished by imposing 
the effective DIM charge to £t the following analytical expression: 

3 

with EjCtj = Zjq — 1. 

As a clear instance of the numerical problems mentioned and the way propose to solve 
them, we show, in Figure [H the orbital (r) of the ground state of the Kr atom (part (a)). 
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Figure 1: (a) Hartree-Fock orbital corresponding to the ground state of the Kr atom. 
It presents two nodes, a genuine one at r 0.06 a.u., and a spurious one at 1.51 a.u. (shown 
in the inset), (b) Dashed line: The corresponding inverted effective charge (r), spoiled 
by the presence of poles. Solid line: Depurated effective charge. 

and its correspondent effective charge (r) (dashed line curve, in part (b)). First, note 
that the 2s orbital has a genuine node at r 0.06 a.u. which produces the hrst pole in the 
effective charge, as shown in the lower graph. The node appears at a relatively low-r value, 
so the corresponding charge (see Eq. (|5])) is not very sensitive to its presence. Therefore, 
it is very easy to eliminate the pole from the effective charge (by just erasing a few points 
around this radius). 

All the bound wavefunctions decay exponentially beyond the last turning point rtp, de- 
hned as the position in which the energy equals the effective potential. At hrst glance, it 
seems that the turning point of U 2 s{r) is located around ^ 0.25, and from that point on, 
the wavefunction should start to decay exponentially. From the numerically point of view, 
r 10 Ttp is a good point to stop the inversion, since beyond there, the effective charge could 
begin to diverge. Thus, one might infer that by erasing the points belonging to the neighbor¬ 
hood of the hrst node, and by stopping the inversion about 10 r^p, the inversion procedure 
will work well. However, the dashed curve in Fig.[T](b) shows a completely unphysical Z^f{r) 
resulting from the inversion. A very careful examination of the (r) orbital function ev¬ 
idences the presence of a spurious node at r = 1.51 a.u., in a region where the amplitude 
of the wavefunction is less than 10“^ times the maximum value (see the inset of Fig. [11(a)). 
Even though this node is completely innocuous for practical matters, it produces devastating 
ehects in the inversion procedure, evidenced by the second huge peak in the Z^f (r) curve 
(see Figure [T](b)). This pole is so big that it ahects a broad vicinity and causes the abrupt 
rising of the effective charge for r > 0.5 a.u.. This is really a surprising result since a priory 
there is no reason to suspect that a negligible oscillation in the tail of the wavefunction would 
produce such a big drawback at small distances. Care must be taken then to discard these 
kind of undesired effects. 
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2.3 Optimization 

The adjustment of the parameters aj and Pj also requires carefull work. The key issue in the 
successful approximation is the region chosen for the fitting: it has to be as large as possible, 
in such a way that Z^™{r) overlaps the inverted Z^^{r) across a broad range, allowing 
an accurate htting procedure, but discarding the points surrounding the nodes. Also, the 
inversion must be halted at a particular (as large as possible) r value, as soon as the amplitude 
of the function is too small. Further on, the inversion procedure diverges. Another issue to 
consider is the self consistency within the computer codes used in the calculations and the 
particular code used to generate the input wavefunctions. To that end, we make sure that 
the same specific numerical grid is used, including the derivatives and integrals at the same 
pivots. The optimization procedure is completed by a number of iteration steps, in which 
the parameters are optimized to give accurate energies and wavefuntions. 

Most density functional approximation methods are based on a variational principle, min¬ 
imizing the density functionals according to energy (others are defined by density). Without 
underestimating its importance, energy is only one of the many parameters that character¬ 
izes a quantum state. Different trial functions (having different forms) can produce, through 
a variational procedure, the same hnal energy. A simple example is given by Bartschat— 
in which two different potentials (one having exchange, the other omitting it) led to produc¬ 
ing very similar and accurate energies of the Rydberg series in several quasi-one electron 
systems. However, a further examination of these potentials shows large discrepancies in 
scattering calculations—. Therefore, in addition to the energy criterion, we have included in 
our optimization method a variational procedure to reproduce accurately the wavefunctions. 
This is achieved by optimizing the mean values (1/r) (which characterize the quality of the 
wavefunction near the origin), and (r) (probing it at longer distances). Furthermore, we 
defined the quantity 

IunFir)u^^{r)dr 

Jp^rir)dr • 

to determine the accuracy of the orbitals generated by the diagonalization of the DIM po¬ 
tentials and the original HF orbitals. 

The effective depurated inversion charge Z^™{r) corresponding to the 2s orbital of the 
Kr atom resulting from the optimization is shown -solid curve- in Figure [T](b). As seen in the 
figure, both boundary conditions are fulfilled (at the origin, Z 2 s 36, and asymptotically 
Z 2 s ^ 1, as stated in Eq. ([6])). 

3 RESULTS 

3.1 DIM Potentials, energies and mean values 

The fitting parameters aj and (dj defining the effective charges Z^™{r) in Eq. ([7]) for the 
noble gases Helium, Neon, Argon and Krypton, are given in Table [T] We have limited the 
aj and fdj to six (about two per shell). For Kr, we would probably need two more since 
there are four shells involved. Having these effective charges, we built the corresponding 
DIM potentials V^™{r). By solving the Schrodinger equation (Eq. ([I])), we obtained the 


6 



Table 1: Fitting parameters for the effective charge Z^™{r) for He, Ne, Ar, and Kr, applying 
Eq. (ED. 


nl 

a 

f3 

nl 

a 


Is 

-0.31745 

5.04372 

Kr Is 

5.49263 

0.884768 


1.31745 

2.50032 


3.94437 

16.8769 


- 

- 


25.5630 

3.10032 

Is 

7.367687 

2.417275 

2s 

9.63120 

0.575832 


1.300360 

0.126396 


1.84650 

25.53280 


0.331953 

13.15820 


23.5223 

4.543350 

2s 

0.297739 

17.99390 

2p 

3.20530 

20.83535 


0.668081 

0.067288 


23.6172 

3.928520 


8.03418 

2.47221 


8.17750 

0.636486 

2p 

1.353049 

8.56948 

3s 

6.52203 

0.547357 


0.335881 

0.464942 


24.4475 

3.657030 


7.311070 

2.090634 


4.03047 

16.61770 

Is 

6.727570 

6.177720 

3p 

23.13135 

4.010523 


4.751090 

1.343560 


3.325360 

20.41890 


5.521340 

0.859981 


8.543290 

0.821218 

2s 

8.90271 

1.09779 

3d 

10.05320 

1.04843 


2.36850 

2.93144 


21.81544 

4.25746 


5.72879 

6.95913 


3.131360 

20.6087 

2p 

4.96956 

6.14455 

4s 

3.65988 

0.49000 


1.48464 

10.86843 


26.4565 

3.17799 


10.5458 

1.30005 


4.88362 

15.2031 

3s 

10.3202 

2.33169 

4p 

7.35713 

1.00142 


4.27115 

7.33678 


24.2321 

3.7309 


2.40865 

0.407463 


3.41077 

22.5680 

3p 

8.43753 

3.49259 





2.18200 

10.8595 





6.38047 

1.07080 
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solutions u^i^{r) and the corresponding energies The comparison between the results 

obtained from the diagonalization of the Hamiltonian with the V^™{r) effective potential 
and the original Hartree-Fock orbitals are presented in Table O It is remarkable that with 
such simple analytical expressions for the potentials we were able to reproduce exactly the 
same energies as the HF method. The only exception is the 4p orbital of Kr, in which 
both calculations agree up to the hfth signihcant hgure. The htting procedure also allows 
to reproduce the original HF wavefunctions with an outstanding degree of accuracy. The 
agreement between the HF orbitals u^f {r) and the solutions can be tested through 

the comparison of the mean values (r) and (l/r), and the computation of quantity 6 dehned 
by Eq. (E]). The mean values agree in about 0.1% while the values of 5 are about 10“^. 

Finally, we calculated the total energy for the ground state of each atom, by using the 
following expression: 


^DIM 


E 





(r) + 




(9) 


where the density = \u^™{r)\‘^. The calculated energies are given in Tabled 

together with the total energies obtained by the Hartree-Fock calculations. The comparison 
shows a notable agreement between both calculations, at about 0.02%. 


3.2 The exchange potential 

Orbital-specihc exchange potentials can be obtained accurately by computing the non-local 
Fock exchange operator. A Erst local approximation can be computed with the average 
exchange charge density proposed by Slater-. Another approximation, proposed by Sharp 
and Horton-, consists in attaining a local potential that approximates the exchange operator 
through a variational procedure that minimizes the energy. There are several other more 
elaborated methods that allow us to obtain local exchange potentials'^—. However, these 
potentials are rather difficult to put in a simple and smooth analytical expression, such as 
Eq.©. 

Due to the fact that the Hartree-Fock method does not take into account the correlations, 
our procedure allowed us to obtain in a rather direct way “exact” local orbital-dependent 
exchange potentials. 


f . ( 10 ) 

r J |r — r I 

where is the total density calculated with the (r) wavefunctions. Figure |2] shows 

the orbital-specihc exchange potentials for the ground states of the four noble 

gases He, Ne, Ar, and Kr, calculated with the depurated inversion method DIM. 

In order to discuss our results, in Fig. [2] we plotted the optimized effective potential 
^OEp(^) developed by Talman- (black dotted lines) for the noble gases. It is well known that 
the OEP method Ends the potential which yields eigenfunctions that minimize the expectation 
value of the Hartree-Fock Hamiltonian. However, although very accurate, it always yields an 
energy above the HF energy. For practical applications the OEP potential works very well for 





Table 2: Total and orbital energies, mean and S values for He, Ne, Ar and Kr atoms obtained 
from DIM effective potentials (upper rows) compared with the original Hartree-Fock values 
(lower rows). 


E 

nl 

e 

(r) 

(1/r) 

(5 

He 

-2.8616 

Is 

-0.917956 

0.927313 

1.687251 

8 X 10-1^^ 


-2.8617 


-0.917956 

0.927273 

1.687282 


Ne 

-128.4978 

Is 

-32.772 447 

0.157491 

9.621450 

2 X 10-6 


-128.5475 


-32.772 443 

0.157631 

9.618 054 




2s 

-1.930 391 

0.891336 

1.640 769 

5 X 10-6 




-1.930 391 

0.892113 

1.632 553 




2p 

-0.850 410 

0.967 755 

1.430 252 

6 X 10-6 




-0.850 410 

0.965 274 

1.435 350 


Ar 

-526.8030 

Is 

-118.610 352 

0.086 015 

17.561606 

2 X 10-6 


-526.8175 


-118.610 350 

0.086 104 

17.553 229 




2s 

-12.322 153 

0.411857 

3.562 264 

2 X 10-6 




-12.322 153 

0.412 280 

3.555 317 




2p 

-9.571466 

0.375 269 

3.449 283 

9 X 10-^ 




-9.571466 

0.375 330 

3.449 989 




3s 

-1.277353 

1.426 944 

0.967005 

9 X 10-6 




-1.277353 

1.422 172 

0.961985 




3p 

-0.591017 

1.668 648 

0.817928 

5 X 10-6 




-0.591017 

1.662 959 

0.814 074 


Kr 

-2752.5365 

Is 

-520.165 467 

0.042 441 

35.483 699 

5 X 10-^ 


-2752.0549 


-520.165 468 

0.042 441 

35.498 152 




2s 

-69.903 081 

0.187181 

7.924 967 

2 X 10-6 




-69.903 082 

0.187256 

7.918 830 




2p 

-63.009 784 

0.161695 

7.874 355 

3 X 10-6 




-63.009 785 

0.161876 

7.868 429 




3s 

-10.849 466 

0.537875 

2.644 610 

2 X 10-6 




-10.849 466 

0.537802 

2.637556 




3p 

-8.331501 

0.542 133 

2.530 080 

2 X 10-6 




-8.331501 

0.542 627 

2.522 775 




3d 

-3.825 234 

0.550 922 

2.276 713 

4 X 10-6 




-3.825 234 

0.550 880 

2.276 940 




4s 

-1.152 935 

1.630 081 

0.808 453 

1 X 10-^ 




-1.152 935 

1.629 391 

0.804188 




Ap 

-0.524186 

1.950193 

0.675 555 

3 X 10-6 




-0.524187 

1.951611 

0.669 219 
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Figure 2: (color online). Orbital-specific exchange potentials and V^gp, for the 

ground state of He, Ne, Ar and Kr. 
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Table 3: Orbital and total exchange energies for He, Ne, Ar, and Kr. 



n 

1 

12 3 4 

Total EAHF 

He 

0 

-0.5129 

-1.0258 -1.026 

Ne 

0 

-3.1106 -0.8620 

-12.1080 -12.11 


1 

-0.6938 


Ar 

0 

-5.8760 -1.9470 -0.5742 

-30.1826 -30.19 


1 

-1.7974 -0.4340 


Kr 

0 

-12.2258 -4.5523 -1.9972 -0.5275 

-93.8525 -93.89 


1 

-4.4305 -1.8401 -0.3906 



2 

-1.5280 



the onter shell. At longer distances, all the nl orbitals have a similar behavior accompanying 
the OEP exchange potential. We noticed that the exchange potentials of the orbitals having 
a common angular momentum I resemble to each other (see Ar for instance). This was 
suggested in a work by Herman et al.— where an /-averaged exchange potential for each set 
of electronic states was calculated as a modihcation of Slater’s average exchange potential. 

According to Eq. flTOj) all orbital-specihc potentials should approach the same value at 
r = 0, since Z^™{r) = TV^™{r) approaches regardless of nl (the second and third term 
are the same for every orbital). However, from Fig. [2] it appears that the potentials for the 
different orbitals approach different values at the origin. This is a consequence of the fact 
that every DIM potential tends to Zjv with different behavior, determined by their htting 
parameters. In fact, for very low r values V^™{r) ^ ~ E‘^(r), but they all have 

strictly the same value at r = 0. 

As a hnal test for our method, we calculated the total exchange energy as given by 


s" = E = 

nl 





( 11 ) 


Table E] displays the orbital exchange energy as well as the total exchange energy for He, Ne, 
Ar and Kr. The total exchange energies are compared with the exact atomic Hartree-Fock 
(EAHF) values given by Becked, with very good agreement. 


3.3 Nitrogen DIM and Exchange Potentials 

The procedure developed here is not limited to noble gases or closed shells. As an example 
we will apply the method to Nitrogen. The lower conhguration 2p^ of Nitrogen gives rise 
to three different terms: 2"^S, 2^D, 2 ^P. Each of them is described by a different electronic 
density. The htting parameters that dehne the term-dependent effective charges are given 
in Table |4] for each of the terms. We built the corresponding DIM potentials from these 
effective charges. By using these potentials we solved the Schrodinger equation (Eq.([T])) 
for every term, obtained the solutions, the energies, and the corresponding mean values (r) 
and (l/r). The comparison between the orbitals obtained from the diagonalization of the 
Hamiltonian with the effective potentials and the original Hartree-Fock orbitals are shown 
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Table 4: Fitting parameters for the effective charge z^r^r) for 2 2 and 2 terms of 

Nitrogen. 


nl 

2 

a 

4g 

/S 

2 

a 

^D 

/S 

2 

a 

"P 

Is 

5.25634 

1.26207 

5.18635 

1.22410 

5.18635 

1.21779 


0.743660 

8.02844 

0.813650 

7.56800 

0.813650 

7.56740 

2s 

2.45281 

3.51271 

0.398100 

0.239738 

0.890660 

0.830615 


0.833570 

3.38654 

1.85412 

1.03105 

3.66999 

3.14946 


2.71362 

0.894699 

3.74778 

2.85313 

1.43935 

0.740427 

2p 

3.64345 

1.24069 

4.01052 

1.28744 

1.89769 

1.16557 


2.05501 

5.35135 

1.85517 

5.70858 

1.77430 

5.68782 


0.301540 

0.286609 

0.134310 

0.267987 

2.32801 

1.40925 


in Table [5l The mean values (r) obtained with the DIM effective potentials agree with the 
HF values in about 0.1%, and the (l/r) mean values agree in about 0.2%. The calculated 
total energies for each term of the Nitrogen atom using Eq. fllip are presented in Table 
O The agreement between the DIM total energies and the original HF total energies is 
excellent, of about 0.04%. Figure [3] shows the nZ-orbital exchange potentials for the 2"^S, 
2 ^D and 2 terms, calculated with the depurated inversion method. Again, to compare our 
results, the exchange potential given by Talman- (OEP) is presented in the hgures in light 
grey. Figure [3]^a) illustrates the exchange potential for the Is orbitals for the different terms, 
showing an overall similarity. The OEP potential behaves like the only at short 

and large distances. Figure EJ^b) shows the exchange potentials for the 2s orbitals. In this 
case, noticeable differences between the term-potentials arise at low values of r. For r higher 
than 0.5 a.u., all the term-potentials become indistinguishable and agree perfectly with the 
OEP potential. A pecularity observed in the hgure is that the OEP potential agrees very well 
with the 2^D term. Figure El^c) displays the 2p exchange potentials, which 

behave similarly for all the terms. However, since the OEP potential is the same for all the 
orbitals and terms, it disagrees completely with the V^pir) at short distances. 

Table [6] presents the total exchange energy and the nZ-exchange energies of the 2 ^S, 2 ^D 
and 2 ^P terms. The Is-exchange energy for all the terms are the same, as expected for a 
closed-shell orbital. Similarly, the 2s-exchange energy varies slightly, with a difference of 
0.08%. However, this is not the case for the 2p-exchange energy, which varies signihcantly, 
having discrepancies of about 18% between the different terms. The total exchange energy 
computed with Eq. flTT]) for the terms are compared with the exact atomic Hartree-Fock 
(EAHF) exchange energy, with an agreement of about 0.1%. 

4 CONCLUSIONS 

A crucial requirement of the density functional method is the accurate representation of the 
exchange functional. On the other hand, the atomic collision community needs accurate one- 
electron potentials in order to generate the bound and continuum states on the same footing 
for further calculations of collisional processes. These potentials need to be worked out for 
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Table 5: Total and orbital energy and mean values for the 2^S, 2^D and 2^P terms of N 
obtained from the DIM effective potentials (upper rows) compared with the Hartree-Fock 
values (lower rows). 


E 

nl 

e 

(r) 

(1/r) 

2^8 

-54.37617 

Is 

-15.629 06 

0.228 30 

6.648 63 


-54.400 93 


-15.629 06 

0.228 30 

6.653 24 



2s 

-0.945 32 

1.33448 

1.080 37 




-0.945 32 

1.332 28 

1.07818 



2p 

-0.56759 

1.412 68 

0.954 98 




-0.56759 

1.409 63 

0.95769 

2^D 

-54.275 57 

Is 

-15.666 39 

0.228 29 

6.649 29 


-54.29617 


-15.666 39 

0.228 26 

6.653 88 



2s 

-0.963 67 

1.32917 

1.086 44 




-0.963 67 

1.326 32 

1.08318 



2p 

-0.508 66 

1.448 78 

0.938 82 




-0.508 66 

1.446 62 

0.942 08 

2^F 

-54.208 56 

Is 

-15.69160 

0.228 24 

6.650 36 


-54.228 10 


-15.69160 

0.228 24 

6.654 30 



2s 

-0.976 34 

1.325 62 

1.08712 




-0.976 34 

1.322 32 

1.086 56 



2p 

-0.471 30 

1.471 76 

0.929 82 




-0.471 30 

1.473 01 

0.93155 


Table 6: Orbital and total exchange energies for 2^S, 2^D and 2^P terms of Nitrogen. 



Is 

2s 

2p 

Total 

EAHF 

2^8 

-2.1175 

-0.4776 

-0.4711 

-6.6034 

-6.596 

2 2d 

-2.1175 

-0.4777 

-0.4262 

-6.4688 


2 2p 

-2.1175 

-0.4780 

-0.3973 

-6.3827 
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Figure 3: (color online). DIM exchange potential for the (a) Is, (b) 2s and (c) 2p 

orbitals, for the 2^S, 2^D and 2^P terms of Nitrogen. 
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any n/-specific orbital, a feature that in general is not present in the chemistry community 
functionals. In the present work we devised and implemented a depurated inversion method, 
which allows to obtain the intended potentials through a very simple analytical expression 
of the effective charges. The method consists in the inversion of a Kohn-Sham equation, 
in which the KS orbitals have been replaced by the Hartree-Fock orbitals. By means of 
diagonalization we have achieved accurate wavefunctions having almost perfect agreement 
with the original Hartree-Fock wave functions. The quality of the potentials obtained by 
the present method is remarkably good. We applied the developed methodology to the 
calculation of the ground state orbitals of noble gases and the Nitrogen atom. It is worth 
mentioning that the same technique can be used for any other level, i.e., it is not limited to 
the ground state. 
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